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Abstract 



We present a new method to derive transport equations for non-relativistic 
quantum many-particle systems. This method uses an equation-of-motion 
technique and is applicable to interacting fermions and (or) bosons in ar- 
bitrary time-dependent external fields. Using a cluster expansion of the r- 
particle density matrices the infinite hierarchy of equations of motion for 
many-particle expectation values is transposed into an equivalent one in terms 
of correlations. This new hierarchy permits a systematic breaking of the hier- 
archy at any order. Diagrams are derived for these transport equations. In a 
second paper the method is tested for exactly soluble electron-phonon models 
in one dimension. 



I. INTRODUCTION 

In this work, the first one in a series of two papers on transport equations, we present a 
new formalism for the time evolution of one-particle distribution functions which allows the 
inclusion of many-body correlations. In the second paper we derive the transport equations 
for an electron-phonon system and compare the numerical results with the exact solutions 
of one- dimensional models (polaron model, Tomonaga-Luttinger model). 

Very often non-equilibrium phenomena in quantum many-particle systems are studied 
by the use of real-time Green's function techniques (Martin and Schwinger 0,0], Kadanoff 
and Baym H, Keldysh |J, Langreth and Wilkins PJ7]]). For a review see Ref. In 
the Keldysh formalism the perturbative expansion uses a time integration over a path in 
the complex plane. If only interested in the one-particle distribution function, one usually 
applies the generalized Kadanoff-Baym ansatz PJ^Q], which was justified by Schoeller [H^ 



as a partial resummation in the diagrammatic expansion. The transport equation obtained 
in the Born approximation is of the form of a non-Markovian Boltzmann equation with 



renormalized one-particle propagators, usually taken the same as in equilibrium. In Ref. flO 
initial correlations are also discussed, but it does not seem to be possible to apply the same 
resummation method including the initial correlations. This stems from the fact that the 
correlations are not taken into account as time- dependent quantities. 
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Motivated by the work of Schoeller, we present a method which enables the inclusion of 
the correlations as dynamical quantities. Only single-time expectation values of operators 
are considered in contrast to Green's function techniques. In the course of this work (cf. part 
II) it will become clear that the transport equations with correlations already include the 
renormalization of one-particle propagators in the dynamics. There is no need to calculate 
retarded and advanced Green's functions independently of the distribution functions. 

As the dynamics of all quantities is just given by the Hamiltonian in the form of an 
equation of motion (EOM) (Schrodinger or von- Neumann equation), it seems to be more 
natural to consider only single-time quantities. The problem to be solved is how to cut 
off the infinite hierarchy of equations of motion in a consistent way. It will be shown that 
the expansion of the expectation values of many-body operators in terms of correlations 
(or, diagrammatically speaking, "connected parts") allows a systematic decoupling of the 
hierarchy at any given order. The lowest order (almost) agrees with the usual decoupling 
procedure to obtain the Born approximation. 



II. EQUATIONS OF MOTION AND CORRELATIONS 

In this section the method will be explained for a fermionic system with two-particle 
interaction in time-dependent external fields described by the Hamiltonian 

# = E + E hf + ± ■ ^ E v hi2 , hj MA^h = H t + v . (i) 

The fermionic annihilation and creation operators are denoted by ipi and ipj, the index i 
referring to the one-particle states. In the interaction term V the matrix elements Vi 1 i 2 j 1 j 2 are 
supposed to be anti-symmetric in the first and second pair of indices (v^ j u - 2 = —v^j^ = 
v hi2,j2ji)- The external fields are described by the time-dependent hermitian matrix h^f"{t). 
Sometimes we use the notation hij(t) = 8ijei + hff 1 ^) for the one-particle Hamiltonian 
matrix. 

The method is not restricted to this type of systems. It may as well be applied to systems 
with bosons, even with "anomalous" (i.e. not conserving the particle number) expectation 
values like (bb). The necessary small alterations are mentioned from time to time in this 
presentation and in part II the method is applied to an electron-phonon system. 

Throughout this work we use the Heisenberg representation, so that the following equa- 
tion of motion for operators holds: 

j t A(t)+i[A,Ht\(t) = -i[A,V](t) . (2) 

For the case of an operator A = ip^ ■ ■ ■ ipk„Wk' • • • tyy and the Hamiltonian of Eq. (HD the 
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equation of motion reads 

(3) 
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* E { h iK (*) (V>fc a • • • ^< • • • 4^1) (*) + • • •} 



E v hi2,jij2 

»1>»2 
Jl J2 



ft 1 J 



or more specifically for the case without external fields 



d 



j t + i (efci + • • • + efc n - e k[ 



(4) 



ii .J2 



The same equations of motion hold for the expectation values of these multi-particle op- 
erators taken with an initial statistical operator p a t time to, (A) t = Tr [p A(t)]. (The 
Heisenberg representation of an operator is supposed to be such that A(t ) = A.) For our 
method to work p does not have to be of any special form. As the case of a non-diagonal 
one-particle matrix hij(t) is notationally cumbersome, but does not pose any additional prob- 
lems in the derivation of the kinetic equations, we will carry out the derivation for the case 
of a diagonal one-particle Hamiltonian, adding the necessary alterations for non-vanishing 
external fields in the end. 

Thus, given the initial state po (mixed or pure) at time to, we have to solve an infinite 
hierarchy of ordinary differential equations (ODE). Of course, this is usually not feasible 
and some kind of decoupling procedure is needed. Starting with one-particle expectation 
values (iplipk')t the EOM involve expectation values of two-particle operators. In the equa- 
tion of these quantities the "collision term" (the right hand side (rhs) of Eq. (p|)) is often 
approximated by factorizing all expectation values into products of one-particle expectation 
values as if the Wick theorem holds for an arbitrary state of the system. In this way a 
non-Markovian Boltzmann equation in Born approximation is obtained. 

In this work we present a method to extend this procedure beyond the Born approxi- 
mation. It yields again an infinite hierarchy of EOM, which is still exact, but now in terms 
of correlations. In contrast to the usual hierarchy of EOM mentioned above, this hierarchy 
allows the decoupling at any given order without any ambiguities. The essential ingredi- 



ent is the following cluster expansion 
annihilation operators, denoted by Bf 



(B 1 ) t = (B 1 y t , 
(B 1 B 2 ) t = (B 1 B 2 y t + (B 1 y t ■ (B 2 y t , 

{BiB 2 B 3 ) t = {B l B 2 B,) c t + {B l B 2 ) c t {B,) c t 



of products of field operators, i.e. creation and 



(5) 

{BrtiB.B^ + {B x )\(B*Btf t + (B 1 ) c t (B 2 ) c t {B 3 r t 



The n-th order correlations (BiB 2 ■ ■ ■ B n )1 are defined recursively. The operators need not 
be normal-ordered. Every expectation value {BiB 2 ■ ■ ■ B n ) t is decomposed into sums over 
products of correlations in the following way: The sum extends over all disjoint partitions of 
the set {!,... ,n}. The product is taken over the correlations corresponding to each subset 
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in the partition. Inside each correlation the order of the field operators is retained. The 
overall sign of the product is determined by the sign of the permutation of fermionic field 
operators necessary to disentangle overlapping correlations. As we assume that the initial 
statistical operator as well as the Hamiltonian conserve the total number of fermions, only 
correlations with an equal number of fermionic creation and annihilation operators occur. 
Because of this property the sign of each decomposition into correlations is well defined. 

So far, the cluster expansion has been written down for both, fermionic and bosonic 
operators. For fermionic systems the particle number must be conserved and the usual rules 
for permutations of fermionic operators apply. The one-particle correlations coincide with 
the corresponding distribution functions 

{^))t={M]Yt > (6) 
<^i>t = <V>i) t c , (7) 

and the two-particle correlations (^i 1 ^i 2 ^j 2 V'j 1 )t; f° r example, are defined by 

^M)^i) t = (M^Ayt + (Mi)$(Mi) c t - (Mir t (Mk) c t ■ (8) 

From the canonical anti-commutation rules it follows immediately that 

(^}) c t = S ij -(^ i ) c t . (9) 
In all higher-order correlations the fermionic operators anti-commute, i.e. 

{...1tf>l--)i = -(-1>]*i-)l , (io) 
(...4> i ij) j ...)t = -(...ij) j ij)i...)t , 

which can be easily proved by induction. Therefore, it is sufficient to consider the following 
anti-normal ordered correlations: 

(^•••^<—^>* ,n>l . (11) 

We will call the one-particle correlations (n = 1) contractions. 

Any expectation value of the form (ip kl ■ ■ -^kjl^y • • •'Wu )t can be expressed in terms of 
correlations as can be seen directly from the definition of the correlations in Eq. (0). This 
enables us to transform the infinite hierarchy of EOM into an equivalent one in terms of 
correlations. So far, nothing seems to be gained. But we will see that the EOM for the 
correlations allow an easy and unambiguous breaking of the hierarchy at any desired order, 
in contrast to the usual EOM. The lowest order is usually of the type of a non-Markovian 
Boltzmann equation in Born approximation. 

In the next section we derive diagrams for the correlations which are as convenient 
as usual Feynman graphs for organizing the calculations. It turns out that the temporal 
evolution of correlations is determined by "connected" diagrams only. 

We start with the EOM for the one-particle distribution function 
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+ i{e kl -e fe /J (tpkM'Xt 



(12) 



n.»2 

Jl ,J2 



*1>*2 
n J2 



+ 



The reader might wonder why in the first step we did not calculate explicitly the commutator 
leaving only terms with four instead of six fermionic operators. It turns out that the cluster 
expansion and the derivation of diagrammatic rules is easier when one leaves the commutator 
as it is. Of course, after the cluster expansion of the commutator the correlations with six 
field operators cancel because of the rule Eq. (|T0|). Only terms with at least one contraction 
remain because the non-trivial part of the fermionic anti-commutation rules shows up only 
in the contractions. In addition, we simplified the result by combining terms differing only 
in the exchange of the two annihilation or creation operators of the interaction term. This is 
possible as we use anti-symmetrized interaction matrix elements. In this way the prefactor 
of 1/4 = 1/2! • 1/2! is sometimes replaced by 1/2! or 1. Later, we will give a simple 
diagrammatic rule to determine this so-called symmetry factor. For pedagogical reasons we 
have not simplified (Vwi)? + { r 4 ) j'4 ) \)t = $ij as we wm do from now on. 

In the next step, the EOM for the newly encountered correlations (V'feiV'fe^I'V'fc')? i s 
needed. It can be calculated in the following way: By definition 



so that the EOM for (4 , k 1 ipk 2 i J kt,' l l J k')t * s nee ded. It reads 
d 



(13) 



dt 



(14) 



h >h 



+ 
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— {k\ <-> k 2 ) — (k[ <-> A; 2 ) + (£4 <-> A; 2 , &i ^ & 2 ) 
+^i,fci(V'l 2 ^ 2 )t(^ii'0Jfc 2 '0fe^i' 1 )t - (h <-> fe) 

-(4^ 1 > t c (4^ 2 ) t c (^ 1 ^> t c (^- 2 ^)0} . 

where we have already used the cluster expansion for the "collision term", i.e. for the rhs 
of the equation. From this equation, Eq. (|T2|), and the definition of the corresponding 
correlation Eq. (|13|) , we obtain the following EOM for (^ki^k^y )* : 



dl 



~ i (^fci^)t 



+ 



dt 



E 

«li»2 
Jl ,J2 



4 (<v> fel ^<«l> t c - (^^(viwf) (^^44) 



(15) 



—(ki «-> k 2 ) — (k[ <-> k' 2 ) + (k\ <-> ^2, fej <-> fc 2 ) 

+ ((^ 1 4)?(V'teV'J 2 )?(4v'i 1 )t(V'i^i 2 )? 
-(4^ 1 ) t c (4^ 2 )?(^ 1 ^)t c (^ 2 ^) t c )} • 

We would like to mention that the diagrams presented in the next section facilitate the 
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determination of the EOM in quite the same way as Feynman graphs do in equilibrium 
theory. Note that the EOM for the correlation is the same as the EOM for the corresponding 
n-point function with just some terms missing. Diagrammatically these terms correspond 
to unconnected graphs. 

Now we want to describe the procedure for the breaking of the hierarchy of EOM. With 
the two-particle interaction the EOM for a given correlation of order n contains correlations 
of order I < n + 1. If a calculation up to order n is desired, we just neglect the correlations 
of order n + 1 in the EOM for the correlations of order n and get a closed system of ODE 
for the correlations of order I < n. This procedure has the advantage of working without 
ambiguities at every given order. In addition, the transport equations are already of the 
form in which they are usually solved, 1.6. cLS db system of ODE. This system allows for the 
inclusion of initial correlations up to order n, that means that additional information about 
the initial state po can be included compared with the usual approaches. It should be clear 
that our approximation method is in some sense an expansion in powers of V ■ At, where V 
is the typical interaction strength and At = t — t is the time passed since the initial time. 
But at the given order it is the best possible approximation for the one-particle distribution 
function. No resummation as in ordinary perturbation theory |10[ is needed. 

The same procedure does not work for the usual hierarchy of EOM because the EOM 
for the 2n-point function (ip^ ■ ■ • ipk n rk' '"Yy)t essentially involves only 2(n + l)-point 
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functions. Therefore, these functions cannot just be neglected. Some kind of decoupling 



mechanism is needed. This was done e.g. by Zimmermann and Wauer [11 in one order 
higher than Born approximation for a Jaynes-Cummings model. But as the only consistent 
way of doing it is in terms of a cluster expansion, it is best to start with correlations from 
the beginning. 

Another noteworthy feature of our method is the fact that the derivatives (4) (V'fcjV'iL )t 



at the initial time t = to are correctly described for < I < n — 1 when the cut-off is at 
order n. A further advantage of using correlations is the fact that they vanish for all orders 
n > 2 in the case of a non-interacting system in a grand-canonical ensemble (Wick's theorem 
fl"2||), in contrast to n-point functions. For this reason the EOM for correlations yield better 
results for systems with weak interaction showing relaxation into thermal equilibrium (cf. 
part II). The most convincing argument in favor of this method are the very promising 
results in comparison with the exact model calculations presented in II for one-dimensional 
electron-phonon systems. In three dimensions the order of the system of ODE is already 
huge at the level of the Born approximation so that it seems not feasible to go even one 
order higher. Even in a homogeneous system with momentum conservation, for example, 
the n-th order correlation has In — 1 free momentum indices, i.e. the number of ODE 
is of the order N^ 2n ~ 1 ' d , where N d is the number of one-electron states. Because of this 
fact numerical calculations in three dimensions are often performed assuming an isotropic 
distribution in fc-space, thus reducing the effective dimension. With electron-phonon instead 
of electron-electron interaction, the situation is slightly better, because the usual electron- 
phonon interaction term contains only products of three instead of four field operators. 

Before turning to diagrams, we discuss the kinetic equations in the lowest approximations. 
Neglecting all correlations of order n>2we obtain 
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j t (^A'/t + ^^)(Mk) c t -<E<^}>*^(*) = ° » ( 16 ) 

3 » 

which is just the propagation in the external fields with the one-particle energies corrected 
by the time-dependent Hartree-Fock energies v^ F (t) := J2i 2 ,j 2 v Hi,3h (V^'fe)*) ^'(0 := % e « + 
h^it) +v^ j F (t). Note that we have added the contribution of the external fields again. This 
equation does not contain a "collision term" . At least the next order has to be retained. We 
obtain the following closed system of ODE: 

^fa^ + iE^W^ ( 17 ) 

3 i 



_ n 

J2,jl,32 11,12,32 



jMk.Ml^lXt +«E{ e *ii( t )^A^A 1 >t c + e fc 2J W(^i^^^)t} (18) 



(ft 



\ji32 n*2 
+ *E ^Wfc^ {t)(Mh$\M0 c t ~ ( fc i ^ ^) - <-> fc 2 ) 



«2J2 



(fci <-> ^2, fc^ <-> A;^ 



2/ 



= -» E ^iiiaJiia ((^ 1 V'l 1 ) t C (V'fe2V'J 2 )t(^ 1 ^ 1 )t(^^2)t 

*1»«2 
31 J2 

- (V'l 1 ^ 1 )t(V'l 2 ^2>*(^ 1 ^ 1 )*(^2^> t C ) • 

The time- dependent Hartree-Fock corrections of the one-particle energies as well as the terms 
with 

v kl k 2 , jlJ2 (t) = E v hi2,3d2 {{MiYM^iYt - (4M c t (^k 2 Y) (is) 

U»2 

and 

W kli2 , k ' lj2 (t) = E v hi2 , jlj2 ({MiM^nYt - (^iMtiMlOt) > ( 2 °) 

which correct the correlated propagation of two particles, are not obtained by a naive de- 
coupling procedure for the two-particle expectation values. 

For the special case of a homogeneous, spin-independent initial state p and a momentum 
and spin conserving Hamiltonian without external fields, the EOM simplify considerably. 
We denote the one-particle distribution function by n k (t) := {ip{ipk)t with k = (k,a). 
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dt 



^ /C2 /C-^ ^2 



— + i(e fcl (*) + e k2 (t) - e k[ (t) - e k ; 2 (t)) 



(21) 



(22) 



3132 



«1«2 

+ ^ 51 V kli2,k[j2 
*2J2 

— (A;x <-> fc 2 ) — (&i <-> fc 2 ) + (fci <-> fc 2 , <-> fc 2 ) 



= -iufcifca,*^ ((1 - "fciWX 1 - n k2 {t))n kli {t)n k > 2 {t) 

- n kl (t)n k2 (t)(l - n^(t))(l - ny 2 (t))) . 

In this form, the ODE are solved numerically. Neglecting the Hartree-Fock contribution and 
the other corrections on the lhs of the second equation, this EOM can easily be integrated: 

(M^kX = (^A^l'X e- i{ ^ +ek ^r%^ (23) 

-i{t kl +€ k2 -e k , -e k , )(t-t') 



-iv klk2 , k > lk ' 2 I dt e 



to 



x ((1 - n kl )(l - n k2 )n k >n k > 2 - n kl n k2 (l - n fc j)(l - ra^)) (if) . 



In the case of no initial correlations, the transport equation for the one-particle distribution 
function can then be written as an integro-differential equation: 



d ft 

—n k (t) = - \ v kk 2 ,k[k' 2 \ 2 / dt' cos((e fc; + e k > 2 - e k - e k2 )(t - t')) 
at k 2 k' 1 k' 2 Jt ° 

x (n k n k2 (l - n fci )(l - n*/) - (1 - n k )(l - n k2 )n k >n k ^ (f) 
For a local two-particle interaction (k\ = (&i,<7i), etc.) 

V klk2 , k ' lk ' 2 = \^ kl - k '^aia'^(T2(j' 2 — ^fc 1 -fc? ( W 2 <W 1 } ^fcj+^^l+fc^ 

and for an initially spin-independent distribution := n^ a (t)) we obtain 



(24) 



(25) 



(26) 



^ ^2 



x (/W 2 (! - "fc-jX 1 " « fc j) - (! - - ^K^Al) (0 > 

which is the usual non-Markovian Boltzmann equation with a collision term for the electron- 
electron interaction in Born approximation. 
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III. DIAGRAMS 



The diagrams described in this section represent the differential equations for the corre- 
lations. Therefore they contain only one interaction vertex in contrast to ordinary Feynman 
graphs. In addition there are graphic elements for the correlations, the so-called correlation 



bubbles ||10|| , which in our context stand for time-dependent quantities and not for initial 
correlations. Note that there is no time or energy integration associated with a diagram as 
it characterizes the temporal evolution of the correlation at time t. 

The diagrams are derived for our example, but they are easily extended for systems with 
bosons and fermions, for systems with three-particle interactions etc. The main elements of 
the diagrams are the following: 

1. The vertex: An internal vertex (Fig. [I]) represents the interaction. The matrix ele- 

ment Vi 1 i 2 j 1 j 2 is associated with it. Incoming (outgoing) lines characterize annihilation 
(creation) operators. In order to simplify the determination of the correct sign, we 
sometimes use another diagrammatic element in which the creation operator ip\ n is 
paired with the annihilation operator ipj n (n = 1,2) (Fig. 0). 

An external vertex (Fig. |3|) represents a single creation or annihilation operator. 

2. The contraction: A contraction (Fig. 0) connects two vertices and refers either to the 

one-particle expectation value (V'ftiVv)* or (— (^fe'^fei)*)- This will be explained later. 

3. The correlation bubble: A correlation bubble possesses incoming and outgoing lines. 

The total number of lines must be different from two, otherwise it is a contraction. 
For the case of a fermionic system only the following correlations occur: 

(^■■■Mk'--4') C t >™> 2 • (27) 

ft i 

The corresponding correlation bubble is shown in Fig. |5|. In order to determine the 
sign, it is better to pair the fermionic creation and annihilation operators in such a 
way that ip^ and %pl, constitute a pair (i = 1, . . . ,n), which is graphically expressed 

i 

as shown in Fig. |5| 

We begin by describing the diagrammatic rules for the part of the time derivative due 
to the interaction of the n-particle expectation values (0^ ■ • -ipkn^l' • • 'Wli)t- The corre- 

n 1 

sponding diagram consists of 2n external vertices and one interaction vertex. These vertices 
are connected in all possible ways by contractions and correlation bubbles. This corresponds 
to the cluster expansion of the commutator ([if>ki ' ' ' ^kJ^X* • • • > V])t- The following rules 
apply for the prefactor and the sign: 

1. The prefactor: Usually the factor 1/2! ■ 1/2! can be omitted due to the symmetry prop- 
erties of the interaction matrix elements v because the graph which is obtained by a 
permutation of the creation and annihilation operators in the interaction term yields 
the same contribution. For this class of graphs only one "unlabeled" graph is written 
down. There is only one exception if the interaction vertex is connected with a cor- 
relation bubble by two lines in the same direction, a so-called pair of equivalent lines. 
In this case a factor of 1/2! remains as the corresponding permutation yields the same 
"labeled" graph. In addition, there is always the factor (— i) of Eq. (0). 
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2. The sign: In addition to the pairing of creation and annihilation operators in the interac- 
tion vertex and in the correlation bubbles, the external vertices are grouped into pairs 
of a creation and an annihilation operator in the same way as for the correlations, i.e. 
ipki and ipy = 1> ■ ■ ■ > n ) f° rm a pair. Then, the overall sign is given by 

i 

{-m-if , (28) 



where L is the number of fermionic loops in the deformed diagrams in which the paired 
external vertices coincide and (—1)™ is the sign of the permutation normal-ordering 
the external pairs. This sign rule can be proved in the same way as for usual Green's 
function methods [|T3]|. 

The description of the resulting differential equation 

d 



dt 

1 



-I 



^ ( 1) + 9 „ e X/ W U«2,jij2^diagr. > (29) 



diagr. n.»2 
31 ,32 



(n e = number of pairs of equivalent lines, L is introduced above), is thus complete except for 
the "collision term" Xdi agr ., which is composed of the contractions and the correlations. In 
the cluster expansion every partition of (ipki ■ • • t/^Vv ' ' 'rk' 'V)t corresponds to a partition 
of (V ■ ipk x ■ ■ '^knYk' ' ' 'rk')t- Since the correlations are (except for the sign) independent 
of the relative order of the creation and annihilation operators, the two terms differ only in 
the contributions of the contractions. Therefore the collision term Xdi ag r. is the product of 
the following contributions: 

1. all correlations; 

2. all contractions which start and end at the interaction vertex in the normal-ordered 
form {—{il , i'4 , j)t) as the interaction term was supposed to be normal-ordered; 

3. the contractions between external vertices in anti-normal order (V'fc-V'I')t as they ap- 
pear in this order on the lhs of the differential equation; 

4. the following contribution of the remaining contractions between the interaction vertex 
and the external vertices, when these lines are labeled as in Fig. [7| 



■ ■ ■ ■ ■ ■ (-(4>^ r )t) (so) 



At least one contraction between the interaction vertex and an external vertex is necessary 
for the graph to yield a non-vanishing contribution. Only in the contractions the non-trivial 
canonical anti-commutation rules show up. 
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In order to make these diagrammatic rules more transparent we give an example. The 
graph shown in Fig. § possesses one pair of equivalent lines. The labeled version shown in 
Fig. p] has two loops if ip ki an d "4>u are paired (i = 1,2). Therefore it gives the following 

i 

contribution to the differential equation for { r ^k 1 4 , k 2 ' l P\' 2 i l J \')t (cf. Eq. (|14D ): 

J'l ,12 

The diagrammatic rules for the correlations are the same as those for the corresponding 
2n-point functions except that only connected diagrams must be considered (linked cluster 
theorem). The fact that the correlations are independent of the order of creation and 
annihilation operators (except for the sign) and that the 2n-point functions do depend on 
this order does not pose any problems for the following reason: Except for the sign, the order 
of creation and annihilation operators is of importance only in contractions connecting two 
external vertices. And these contractions do not appear in connected diagrams. 

The linked cluster theorem can be proved in the following way. To each graph for a 2n- 
point function belongs a partition of the external vertices into connected components, e.g. 
the graph shown in Fig. [L0| gives the partition {^ k ^wl> 3 ^k 2 ^k 3 }, {^l^^kA- 

The diagrammatic partition determines exactly one term in the cluster expansion of the 
2n-point function. In addition, there is one and only one connected component containing 
the interaction vertex, in the example it is {fyy , ipki}- We call this component the marked 
component. The term in the cluster expansion with this marked component determines one 
term in the temporal evolution of the 2n-point function. In order to reduce the formalism 
this is explained for the example: 



d , 

+ i [e kl H h e/c 4 - e fc / e k 



dt 



{Tp kl ---ipkA>--^k)t (32) 



+ 



d m ■( 



i.e. the time derivative and the one-particle Hamiltonian act on the marked component 
according to Leibniz's rule. 

We will prove the linked cluster theorem by induction on the order of the correlations, 
corresponding to their recursive definition in Eq. @. Obviously only connected diagrams 
contribute to the time evolution of the contractions as at least one external vertex must 
be contracted with the interaction vertex (and the other one must be connected to the 
interaction vertex because of the fermionic particle conservation). We consider now a cor- 
relation (ip kl ■ ■ -tpkniph ■ ■ 'Wk>)t °f or der n > 2 and assume by induction that the linked 
cluster theorem holds for all correlations of order less than n. The set of all diagrams for 
the corresponding 2n-point function (ip kl ■ ■ -ipkn^h ' ' 'rl')t is decomposed into classes of 

n i 

diagrams with the same partition of the external vertices and the same marked component. 
To each of the classes belongs a term in the cluster expansion of the differential equation, as 
described above. By induction, this term is described by the corresponding class of diagrams 
if the marked component is of order less than n. But the total of these diagrams are all 
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unconnected diagrams. Thus, the remaining term in the cluster expansion of the differential 
equation 



d 



ejfei + l-ej^- £k[ 



(fai ■ ■ ■ ^kn^L ■ ■ ■ V4 Yt 



(33) 



is represented by the connected diagrams. This completes the proof. Of course, it is neces- 
sary to check the sign and the prefactor as well. As far as the factor (—i)l/2 ne is concerned, 
this is clear because the factor belongs to the component containing the interaction vertex. 
We will not give a formalized proof for the correctness of the sign, which should nonetheless 
be apparent as the sign only reflects the necessary permutation to bring the fermionic oper- 
ators into the given order. We close with the remark that the proof is the same for the case 
of a non-diagonal one-particle Hamiltonian foy(t) and that for bosons the induction starts 
with the correlations of lowest order (13"')% and no argument about particle-conservation is 
needed. 

The fact that only connected diagrams contribute to the time evolution of correlations is 
explained by the definition of the correlations because we allowed for all possible partitions 
of the creation and annihilation operators. Therefore, it is necessary to define the cluster 
expansion for bosonic systems with non-vanishing expectation values (B^) t as in Eq. ([|) 
and not as sometimes defined: 



(Bi)t 
(B 1 B 2 ) t 

(B 1 B 2 B 3 ) t 



(BxY t ■ 

(BiB 2 y t 
(BiB 2 y t 



(34) 



[B 3 ) c t + (B 2 )t(B 1 B 3 )l + (B 1 ) c t (B 2 B 3 )t + (B 1 B 2 B 3 )1 



Otherwise, also unconnected diagrams give contributions to correlations. 

We now return to the example given in section [TI| and present the diagrams for the 
correlations of one, two, and three particles. The diagrams for the one-particle distribution 
function representing the Hartree-Fock contribution and the "real" collision term in Eq. flT2p 
are shown in Fig. 
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In order to greatly reduce the number of diagrams, we draw the 
diagrams for higher-order correlations unlabeled and without arrows. These diagrams are to 
be seen as representatives for all the diagrams which arise by labeling the external vertices 
in all possible ways and attaching arrows consistently with the conservation of the particle 
number. The last two graphs in Fig. O, for example, are represented by the single graph 
in Fig. 0. With this simplification there are four types of diagrams for the two-particle 
correlations (cf. Eq. ([15])) shown in Fig. |13[ The numbers of resulting diagrams with external 
labels and arrows are indicated in parentheses. For the three-particle correlations we show 
the diagrams in Fig. H without writing down the corresponding equations. To cut off the 
hierarchy at this order the four-particle correlations are ignored as described in section ||. 



IV. SUMMARY 

We have described the general formalism for a new method of deriving transport equa- 
tions for quantum many-body systems. It is based on the equation-of-motion technique and 
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a cluster expansion of many-particle expectation values. Combining these two elements one 
obtains an exact infinite hierarchy of EOM in terms of correlations which allows a consis- 
tent, unambiguous breaking at any order in contrast to the usual hierarchy. The transport 
equations consist of a set of ordinary differential equations for the one-particle distribution 
functions and many-particle correlations up to a given order. By the initial values of these 
dynamical quantities information (initial correlations) of the initial statistical operator en- 
ters the quantum kinetics of the system. In contrast to real-time Green's function methods 
our method involves only single-time quantities from the beginning. There is no need for 
a generalized Kadanoff-Baym ansatz and the separate determination of retarded/advanced 
Green's functions. A diagrammatic representation of the equations of motion for the corre- 
lations was derived. It simplifies the calculations in the same way as Feynman diagrams do 
in equilibrium theory. The definition of the correlations by a cluster expansion leads to a 
linked cluster theorem for the corresponding diagrams. 

In this first part of a series of two papers we developed the general formalism. The 
application of the transport equations to physical systems is left to the second part. 
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FIGURES 




FIG. 1. An interaction vertex. 




FIG. 2. An interaction vertex with paired creation and annihilation operators. 
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FIG. 3. External vertices. 
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FIG. 4. A contraction. 




FIG. 5. A correlation bubble representing (ip^ ■ ■ ■ tpknVk' ■ ■ ■ )t- 

n \ 
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FIG. 6. A correlation bubble with paired creation and annihilation operators. 
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to corr. 

FIG. 7. The graph for the explanation of the Pauli factors for a general n-body vertex. There 
are r + s lines connecting the interaction vertex with external vertices. The remaining lines go to 
correlation bubbles. 




FIG. 8. A graph contributing to the 4-point function. 
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FIG. 9. A labeled version of the graph in Fig. |8| with paired creation and annihilation operators. 





FIG. 10. An unconnected graph for the 8-point function. 
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FIG. 14. The diagrams for the three-particle correlations. 
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